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We discuss the PPN Solar-System constraints and the GW stochastic background considering some recently 
proposed f(R) gravity models which satisfy both cosmological and stability conditions. Using the defini- 
tion of PPN-parameters 7 and /3 in terms of /(7?)-models and the definition of scalar GWs, we compare 
and discuss if it is possible to search for parameter ranges of /(i?)-models working at Solar System and 
jJll ' GW stochastic background scale. 
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! 1 Field equations and viable f(R) -model 

Let us start from the following action (see [1]) 

1 

s = s g + s m = -j-j / d 4 x^d [R + f(R) + Cm] , (1) 

where we have considered the gravitational and matter contributions and k 2 = 167rG. The non-linear 
f(R) term has been put in evidence with respect to the standard Hilbert-Einstein term R and C m is the 



perfect-fluid matter Lagrangian. The field equations are 

^ : \g^F(K) - R^F'(R) - g^DF'(R) + V^ V F'(R) = ~T$\ 



(2) 



2 » 

Here F(R) = R+f(R) and is the matter energy - momentum tensor. Action (fTJ can be recast in a 
scalar-tensor form. By using the conformal scale transformation g M „ — > e° with a = — In (1 + f'(R)), 
the action can be written in the Einstein frame as follows [2]: 

Se = ^J dAx ^~9 (R - \a pa ^ P ^aO - V(a)j , (3) 

where 



W = e - S («-)-«'/( S («-)) = i-||. (4) 
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The form of g (e~ cr ) is given by solving a = — In (1 + f'(R)) = InF'(R) as R = g (e~ CT ). The 
transformation g M „ — > e a g^ v induces a coupling of the scalar field a with matter. 

Let us consider now a class of f(R) models which do not contain cosmological constant and are explic- 
itly designed to satisfy cosmological and Solar-System constraints in given limits of the parameter space. 
In practice, we choose a class of functional forms of f(R) capable of matching, in principle, observational 
data. Firstly, the cosmological model should reproduce the CMBR constraints in the high-redshift regime 
(which agree with the presence of an effective cosmological constant). Secondly, it should give rise to an 
accelerated expansion, at low redshift, according to the ACDM model. Thirdly, there should be sufficient 
degrees of freedom in the parameterization to encompass low redshift phenomena (e.g. the large scale 
structure) according to the observations. Finally, small deviations from GR should be consistent with Solar 
System tests. All these requirements suggest that we can assume the limits 

lim f(R) = constant, lim f(R) = 0, (5) 

R^oo R^O 

which are satisfied by a general class of broken power law models, proposed in [3], which are 



. In 

R 



F(R) =R- XR C - v c / n (6) 



(*) 



1 



where parameters {n, A, R c } are constants which should be determined by experimental bounds. 



2 Constraining f(R) -models by PPN parameters 



The above model can be constrained at Solar System level by considering the PPN formalism. This ap- 
proach is extremely important in order to test gravitational theories and to compare them with GR. As it is 
shown in [4,5], one can derive the PPN -parameters 7 and (3 in terms of a generic analytic function F(R) 
and its derivative 

A 4 



7 



1 



F'(R) ■ F"(R) 



2F'(R) + 3F"(R) 2 



dR 



(7) 



F'(R) + 2F"(R) 2 ' 

These quantities have to fulfill the constraints coming from the Solar System experimental tests summa- 
rized in Table I. They are the perihelion shift of Mercury, the Lunar Laser Ranging, the upper limits coming 
from the Very Long Baseline Interferometry (VLBI) and the results obtained from the Cassini spacecraft 
mission in the delay of the radio waves transmission near the Solar conjunction. 



Mercury perihelion Shift 


|27 -/3- 1| < 3 x 10" 3 


Lunar Laser Ranging 


4/3-7-3= (0.7 ±1) x 1CT ;5 


Very Long Baseline Interferometer 


|7- 1| < 4 x 10- 4 


Cassini Spacecraft 


7-1= (2.1 ±2.3) x 1(T 5 



Table 1 Solar System experimental constraints on the PPN parameters. 



By integrating last equations (0, one obtains f(R) solutions depending on 7 and (3 which has to be 
confronted with ^ exv and f3 exp plug into such equations the model and the experimental values of PPN 
parameters and then we will obtain algebraic constraints for the phenomenological parameters n and A 
. Determining the value from the equation for de Sitter solutions according to the stability conditions 
F'(R) > and F"(R) > 0. Finally we obtain a good sets of parameters for the model [2]: 

• = 3.38 , n= 1 , A = 2 

• € = V3,n = 2,A>^ 



Copyright line will be provided by the publisher 



adp header will be provided by the publisher 



3 



^(/) = -=, Pc = 3 (8) 



3 Constraining f(R) -models by stochastic backgrounds of gravitational 

waves 

Also the stochastic background of GWs can be taken into account in order to constrain models. This ap- 
proach could reveal very interesting because production of primordial GWs could be a robust prediction 
for any model attempting to describe the cosmological evolution at primordial epochs. The main charac- 
teristics of the gravitational backgrounds produced by cosmological sources depend both on the emission 
properties of each single source and on the source rate evolution with redshift. It is therefore interesting to 
compare and contrast the probing power of these classes of /(i?) -models at hight, intermediate and zero 
redshift [6,7]. To this purpose, let us take into account the primordial physical process which gave rise to a 
characteristic spectrum fl sgw for the early stochastic background of relic scalar GWs by which we can re- 
cast the further degrees of freedom coming from fourth-order gravity. This approach can greatly contribute 
to constrain viable cosmological models. The stochastic background of scalar GWs can be described in 
terms of a scalar field $ and characterized by a dimensionless spectrum We can write the energy density 
of scalar GWs in terms of the closure energy density of GWs per logarithmic frequency interval as ( [8]) 

1 dp sgw _ 3Hq 

y c d\nf ' Pc ^8^G 

is the critical energy density of the Universe, Ho the today observed Hubble expansion rate, and dp sgw 
is the energy density of the gravitational radiation scalar part contained in the frequency range from / to 
/ + df. We are considering now standard units. 

The calculation for a simple inflationary model can be performed assuming that the early Universe is 
described by an inflationary de Sitter phase emerging in the radiation dominated era [9]. The conformal 
metric element is 

ds 2 = a 2 (r])[-dT] 2 + die 2 + h lilJ {7 1 ,^)dx^dx 1 '}, (9) 
and a GW with tensor and scalar modes in the z+ direction is given by [10] 

V,(t -z) = A+(t- *)e« + A x (t - z)e$ + *(i - 2)e« . (10) 
The pure scalar component is then 

v = Ht } > en) 

(s) 

where e^i is the polarization tensor. At lower frequencies, the spectrum is given by 

n iBW (/)oc/- a . (12) 

It is interesting to calculate the corresponding strain, where interferometers like VIRGO, LIGO and LISA 
reach a maximum in sensitivity. The well known equation for the characteristic amplitude [8], adapted to 
the scalar component of GWs, can be used. It is 

,_ib ( 



$ c (/) ~ 1.26 x 10- 18 \-j-J y/hl 00 Cl sgw (f), (13) 

and then we obtain the values in the Table [2] 

At this point, using the upper bounds in Table|2] calculated for the characteristic amplitude of GW scalar 
component, let us test the /(i?)-gravity models, considered in the previous sections, to see whether they 
are compatible both with the Solar System and GW stochastic background. 

Before starting with the analysis, taking into account the above discussion, we have that the GW scalar 
component is derived considering 

$ = --, ( 7 = -m(l + /'( J R))=lnF'( J R), 8a = fff ifl. (14) 

Then we obtain a good sets of parameters for the model [2]: 
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$ c (100Hz) < 2 x 1CT 2B 


LIGO 


<P c {100Hz) < 2 x 1CT 25 


VIRGO 


$ c (100Hz) < 2 x 1CT 21 


LISA 



Table 2 Upper limits on the expected amplitude for the GW scalar component by ground-based-interferometers 
LIGO- VIRGO and space-interferometer LISA. 

• JL = 3.38 , n = 1 , A = 2 

• € = V3,n = 2,A>^ 

such sets of parameters are same as bounds coming from the PPN and some sets reproduce quite well both 
the PPN upper limits and the constraints on the scalar component amplitude of GWs. The results indicate 
that self-consistent models could be achieved comparing experimental data at very different scales without 
extrapolating results obtained only at a given scale [2]. 

4 Conclusions 

The interesting feature, and the main result of this paper, is that such sets of parameters are not in conflict 
with bounds coming from the cosmological stochastic background of GWs. In particular, some sets of 
parameters reproduce quite well both the PPN upper limits and the constraints on the scalar component 
amplitude of GWs. 

Far to be definitive, these preliminary results indicate that self-consistent models could be achieved 
comparing experimental data at very different scales without extrapolating results obtained only at a given 
scale. 
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